Tunable multiple layered Dirac cones in optical lattices 
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We show that multiple layered Dirac cones can emerge in the band structure of properly addressed multicom- 
ponent cold fermionic gases in optical lattices. The layered Dirac cones contain multiple copies of massless 
spin- 1/2 Dirac fermions at the same location in momentum space, whose different Fermi velocity can be tuned 
at will. On-site microwave Raman transitions can further be used to mix the different Dirac species, resulting in 
either splitting of or preserving the Dirac point (depending on the symmetry of the on-site term). The tunability 
of the multiple layered Dirac cones allows to simulate a number of fundamental phenomena in modern physics, 
such as neutrino oscillations and exotic particle dispersions with E ~ p N for arbitrary integer N. 



PACS numbers: 37.10.Jk, 71.10.Fd, 71.15.Rf, 14.60.Pq 

Introduction. The concept of a quantum emulator was first 
introduced by Feynman as a way to avoid the difficulty of 
simulating quantum phenomena with classical computers [1]. 
The idea was to use one controllable system to simulate an- 
other, possibly computationally intractable system. Nowadays 
Feynman's intuition is being implemented in various setups 
and among them, cold gases of neutral atoms play a central 
role [2-4] . Systems of trapped ultracold atoms in optical lat- 
tices have proven to be a remarkable tool for simulating a vast 
range of condensed matter, and lately also high energy physics 
phenomena. Graphene and topological insulators have re- 
cently attracted a significant attention due to their massless 
low-energy excitations, termed Dirac fermions, which govern 
the electron transport properties [5, 6]. They are examples of 
quasi-relativistic dynamics in a non-relativistic environment. 
In an effort to extend such behavior to quasi-particles with ar- 
bitrary spin, a quantum simulator of the so called Dirac-Weyl 
fermions with arbitrary large spin [7] has been proposed. This 
can be implemented by using multicomponent cold fermionic 
atoms trapped in an optical lattice, and by tuning the hopping 
matrices between these internal states according to su(2) al- 
gebra representations. This setup allows us to assign an ar- 
bitrary spin to the emergent low-energy excitations, whose 
band structures display multiple layered Dirac cones at four 
isolated points in the first Brillouin zone, termed Dirac points. 
Such band structure gives rise to exotic properties, e.g., a rich 
anomalous quantum Hall effect and Klein multi-refringent 
tunnelling. The shape of the cones, i.e., the different Fermi 
velocities of quasi-particles are completely fixed by the su(2) 
representation. In this Letter, we show how to relax the above 
constraint, providing a playground for simulations of high en- 
ergy phenomena such as neutrino oscillations [8] in Lorentz 
and CPT breaking/no breaking scenarios [9-11], and modifi- 
cations of dispersion relations [12, 13]. 

The model. For an optical square superlattice filled with 
multicomponent alkali fermions, the Hamiltonian is [14] 

H = ^+^ = -^^? v [T v ] T , T 4 +V;T ,c rT +H.c. 

+LLi°}^ci Tl c rT , (i) 



where H t describes the nearest-neighbor hopping dynamics, 
H the on-site dynamics, and cl T (c rT ) are fermionic creation 
(annihilation) operators with an internal index T. The optical 
potential is assumed to be deep enough so that free hopping 
due to kinetic energy is suppressed, and all system dynamics 
is induced by Raman transitions. In particular the hopping 
along the V = x, y direction is described by a spin-dependent 
operator T v while the on-site spin-flipping is described by O. 
For a detailed discussion of experimental realizations we re- 
fer the reader to Refs. [7, 14]. In such laser-assisted hopping 
schemes, each of the matrix elements of the hopping operators 
is realized via an effective four-photon process, where the spin 
of the atom can be flipped with at most | Amp | = 4 [14]. As ex- 
plained below, we demand non-zero elements of the hopping 
matrices only around the diagonal, i.e., with a superdiagonal 
and a subdiagonal which involve only \Amp\ = 1. 

Tunneling dynamics and multiple layered Dirac cones. 
In momentum space H t reduces to H = ^^W/ZW^W 
where H(k) = — £ v 2r v T v cos(& v ), the spinor *F(Jfc) = 
(eft, ...,C£ n ) T contains the fermionic operators, and k y 
is in units of the lattice spacing. When the hopping 
operators T v are tuned according to the ^-dimensional 
representation of the su(2) Lie algebra, namely, T x = E x 
and T y = E^ which fulfill the corresponding algebra 
[E Z ,E±] = ±2E± and [E^Ey] = 2iL z where E± = E x ± iL y , 
then the low energy excitations of the system are Dirac- 
Weyl fermions with integer or half-integer spin depending 
on whether n is odd or even. A simple representation of 
E x and E^ for spin s with n = 2s + 1 which generalizes 
the Pauli representation for spin 1/2, can be expressed 
in terms of a (n — l)-vector p n = (pi,p2,--- ,p n -i) with 
Pj = y/j(n-j) such that E x = {superdiag{p 7 },subdiag{p 7 }} 
and E y = {superdiag{— zp 7 }, subdiag{/p 7 }}, j = 1, ... ,n — 1. 

It is worth noticing that for any real vector 
p, T x = {superdiagjpjj^subdiaglPj}} and = 
{superdiag{— zp 7 },subdiag{zp 7 }}, the location of the 
Dirac points of H(k) does not change because the spectrum is 
E(k) = £p\gk\ where £ p are the eigenvalues of the T v matrix 
[15], and g k = (2t x cosk X: 2t y cosk y ,0). This observation is 
crucial as the leitmotif of the Letter is to preserve the stability 



of the Dirac points while relaxing the integer or half-integer 
spin structure such that we can construct arbitrary effective 
speeds of light in the emergent quasi-relativistic scenario. 
With this extension, properties that depend on the topological 
aspects of the system such as the topologically invariant 
Berry phase, the topological charge and the Hall plateaus 
will consequently not change, while properties depending 
on local aspects of the system such as the butterfly nature of 
the spectrum and the Klein multirefringent tunnelling will 
change [7] . The resulting spectrum at each Dirac point will be 
a collection of Dirac fermions with tunable effective speeds 
of light and a topological charge N equal to the number of the 
layers, i.e. N = s + 1/2 for half-integer spin. 

As an example, we consider a double layered Dirac cone 
structure with tunable Fermi velocity for each component. We 
parametrize p 4 as p(sin0cos(p,sin0sin(p,cos0). The corre- 
sponding spectrum is of the form £ = p y/i ±x/2/ y/2, where 
X = y3 -f cos2<p + cos40 — cos2(pcos40. Interestingly, a 
birefingent breakup of the doubly degenerate Dirac cones into 
cones with different speeds of light have also been discussed 
in quite a different setup with cold atoms [16]. It is worth 
stressing that our setup allows a breakup, or indeed a coali- 
tion, of any number of Dirac fermions. For example, when 
X = the two layered cones collapse into a single degenerate 
one with topological charge N = 2. This is also the mecha- 
nism for mixing the different Dirac species when on-site dy- 
namics is introduced. The N = 3 situation, on the other hand, 
can be used to mimic the three families of fermions in particle 
physics. These hopping matrices considered above can also 
be seen as generalized spin-orbit coupled systems (see [17] 
for a recent realization in a Bose-Einstein condensate). 

On-site dynamics and Topological Phase Transitions 
(TPT). While the hopping Hamiltonian H t in Eq. (1) allows us 
to create a collection of Dirac fermions at the same location 
in momentum space with tunable Fermi velocities, a constant 
on-site term H Q gives an effective mixing of the different Dirac 
species. This mixing mechanism is fundamentally different 
compared to the mixing mechanism in graphene. In graphene 
the Dirac fermions are located at different sites K± in the mo- 
mentum space. A commensurate perturbation G = K+ — K- 
introduced by lattice distortions is needed to mix the Dirac 
fermions (see [18, 19] for a discussion of the chiral mixing 
in graphene by Kekule texture when constructing the chi- 
ral gauge theory of graphene, and [20, 21] for arbitrary spin 
quasi-particles). Furthermore, opening a gap in monolayer 
graphene is difficult. Our setup allows mixing of any number 
of Dirac fermions of the same chirality, and to open up gaps 
simply by the on-site microwave Raman transitions. 

We now consider the manipulation of triple degenerate 
Dirac cones with a topological charge of N = 3 at a sin- 
gle Dirac point K = |(1,1), i.e., at the same location in 
momentum space. We write the constant on-site term - a 
special case of the O of Eq. (1) - in the form of h l • G 
where h l = (h l x ,h l yJ h l z ), i = 1,2,3, are constant vectors, and 
G = (g x , c-y, G z ) is a vector of Pauli matrices. The unperturbed 




Figure 1: Setup and typical band structures considered in the 
text, (a) Schematic illustration of the optical lattice and its tunneling 
mechanism, (b) The triple layered Dirac cones with three different 
effective speeds of light, (c) The splitting of triple degenerate Dirac 
cones, (d) Gap opening of triple degenerate Dirac cones. The band 
structures in (c) and (d) are used to simulate the exotic and traditional 
neutrino oscillations as discussed in the text. 



case h l = corresponds to p 6 = (1,0, 1,0, 1). The result is a 
combination of mixing and splitting of the Dirac points, which 
are the paradigm of topological quantum phase transitions. 
The corresponding Hamiltonian near the Dirac point is 



H m (k) = 



h u 
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The spectrum is given by Ei = ±\g k — h l \. This expression 
provides a convenient way to manipulate each of the three 
Dirac species. While h l x and h l y control the positions of the 
split Dirac points, h\ controls the gap of the energy spectrum 
(see Fig. 1). By requiring E = and setting h\ = for the 
time being, the split Dirac points are determined by the condi- 



tions 2/rCOS&r — hi 



and 2t y cosky 



- h l = 0. In principle, h l 



may be tuned to appropriately large values to create marginal 
Dirac points with topological charge N = by merging the 
Dirac points of opposite chirality. The applications of such a 
scenario go beyond the realm of condensed matter [22] . How- 
ever, in what follows we will focus on perturbative splitting. 

Neutrino oscillations. Neutrino oscillations (NO) are 
considered by many as a possible window for physics beyond 
the Standard Model. The most accepted mechanism to 
explain such oscillations is in terms of the nontrivial mass 
term matrix. The Lagrangian density describing the flavour 
states vj = (v e ,v^,v T ) with a mixing mass term M is 
jgf = Vf(x)(if i d^ l -M)Vf(x). The flavour of the neutrino 
oscillates as it propagates since the flavour eigenstates do not 
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Figure 2: Anisotropic quasi-neutrino oscillations and T- violations in 
optical lattice with mixing angles (0i2, #13, #23) = (1/4, 1/4, 1/4) ;r 
and splitting h = 0.01 (2tf) (1,0,0). (a) The oscillations of the prob- 
ability P(v e — >• Vfi) with the directional vector p of the momen- 
tum where the CP-violation phase is 8 = . (b) The T-violation 
behaviour with S = n/2 (red), 3n/2 (blue), 0,7T (green), where 
Ar = (P(y e -+ v M ) -P(v M v e ))/(P(v e v M ) +P(v M v,)). 



coincide with their mass eigenstates with a definite mass and 
energy. In general, the weak interaction flavour eigenstates 
can be represented as a coherent linear superposition of 
the mass eigenstates, V/ = Uv m where v T m = (Vi,V2,V3) 
describes the state with definite masses and U is the unitary 
transform matrix that diagonalizes the mixing mass ma- 
trix M, known as the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) lepton mixing matrix. For all currently observed 
NO, the corresponding masses are less than leV and the 
energies are at least E > lMeV, with a Lorentz factor greater 
than 10 6 in all cases. In this ultrarelativistic limit the en- 
ergy is given by E[ = (p 2 c 2 + mjc 4 ) 1 / 2 ~ E + mjc 4 /2E, 
E = c\p\. The time-dependent mass states are 
|Vf(0> = exp[-iEit/h]\Vi), thus v f (t) = UTU*V f (0) 
where T = &a.g{exp[—iEit/h],exp[—iE2t/h],exp[—iE3t/h]} 
is the time evolution matrix, i.e., P(v a —> v^) = 

m 2 c 4 

\Y,jU a jUpjCxp[—i^^-t]\ 2 , which shows that the neutrino 
flavour changes with time due to the mass difference between 
the mass states [8]. The question of which mechanism 
generates such mass is still not completely settled. 

NO as TPT in optical lattices. An analogue of NO can 
be engineered in optical lattices (for a classical analogue of 



neutral-meson oscillation see [23]). The simple model of Eq. 
2 is the starting point to capture the essence of 3 species' 
mixing. Interestingly, our setup allows us to simulate both 
traditional NO and exotic anisotropic NO. In the following, 
we will consider an ultra-relativistic limit for coupled pseudo- 
particles, i.e., \h\ <C \gt\ ~ Q\p\, where p = (K — k). For the 
traditional NO, the analogue of a mass term can be reproduced 
by setting h l = (0, 0, h[) for which the spectrum E { = \g k — h l | 

reduces to E\ ~ c\\p\ + (h[) /(2c/|p|), thus h[ represents the 
mass of the neutrinos. Moreover, our setup also allows us 
to simulate exotic anisotropic NO where a direction depen- 
dent mass can be achieved in the limit of small h l = (h l x ,h l y: 0) 
which we shall focus on in this study. In this case, the spec- 
trum of E[ — \g k — h l \ reduces to E{ — ci(\p\ + p - h l ), where 
p = p/(ci\p\). For sake of simplicity, we consider h l = 
and h 2 = h = —h 3 . The H m (k) in Eq. 2 therefore plays 
the role of the block diagonal Hamiltonian in the mass ba- 
sis. In fact, to observe NO in the lab we have to imple- 
ment the Hamiltonian Hf(k) governing the flavour pseudo- 
particles. In terms of the PMNS mixing matrix U, Hf(k) = 
(t/ 1 ®I 2 )H m (k)(U <8>/ 2 ) = h ®gk -G+M®h-G, where I n 
is the n x ^-identity matrix, and M = t/^diag{0, — 1, l}t/, i.e., 
Mij = UyU^f — UijU^i- In Lagrangian terms, our 2+1 model 
is jSf = ¥f (x)(iy ti dn -MY i h {1 ) x ¥f(x) near the Dirac point. 
Thus the flavour Hamiltonian fits in the family of Eq. 1, 
Hf (k) = H t (k) + H (h) , with H t (k) = h®g k -a and H Q (h) = 
M&h - G. Remarkably, while H t (k) can be implemented by 
the triple-layered degenerate Dirac cones as discussed above, 
the PMNS mixing matrix is completely encoded by the on- 
site Hamiltonian H (h), which allows us in principle to simu- 
late any mixing angles and CP- violating phase - the standard 
parametrization of the PMNS can be found for instance in [8] 
- by the on-site microwave Raman transitions [14]. As an ex- 
plicit example, we discuss the simplified case when all the 
mixing angles are equal to 71 / 4 with 
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and 5 is the CP-violating phase. The quasi-neutrino oscilla- 
tions are shown in Fig. 2. In Fig. 2 (a), the oscillation prob- 
abilities for 5 = are plotted against p, showing anisotropic 
and energy-independent behaviors. 

The Hamiltonian of Eq. 3 is real when 5=0, hence invari- 
ant under time-reverse symmetry T. Since the oscillation pe- 
riod is inversely proportional to the splitting, one obtains for a 
small splitting parameter \h\ ~ 0.01(27r/a), a period which is 
of the size of typical lattices. In high energy physics, evidence 
for NO is presented over a distance L > 100km and no evi- 
dence for oscillations for L < Ikm [8]. In Fig. 2 (b), the effect 
of 5 7^ is considered. The Hamiltonian becomes complex 
and T is violated. Due to CPT invariance, T-violation is equiv- 
alent to CP-violation. The Time-violation behaviour of the 
transition probability P^a -> ^p) ^ P( x ^ f f5 -> ^a) is shown 
for 5 = {0,1/2,1,3/2}^. The signature of CP-violation is 
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rather spectacular. In the optical lattice experiment, direct evi- 
dences of the above phenomena can be obtained by measuring 
the different populations in different points of the lattice by 
the colour resolution or individual atom detection techniques 
[24, 25]. Lorentz and CPT conservation hold if h is covariant. 

Lorentz and CPT violation in optical lattices. Different 
sources of Lorentz and CPT violating oscillation terms can 
be considered in optical lattices. The easiest is to simulate 
anisotropic Minkowski space-time, i.e., t x ^t y , thus c x ^ c y . 
Alternatively, any isotropic and anisotropic Lorentz violating 
mass term [9], J£ = m + im^ + + ^75/^ -f i^ v <5 MV , 
can be conveniently engineered. In the former case m and m$ 
are Lorentz violating if they have an energy dependence [26], 
and their CPT behavior depends on their origin. Phenomeno- 
logical models including terms such as m °c E p , with p > 0, 
are very appealing, and induce strong modifications of the os- 
cillations at high energy. In optical lattices, such terms can be 
induced by time-dependent Bragg scattering. A detail study is 
left for a future work. In the latter case, for constant coupling, 
H is CPT conserving but Lorentz violating whereas the a and 
b terms are both CPT and Lorentz violating. In principle all 
of them can be reproduced by the more complicated on-site 
spin-flipping Hamiltonian as in Eq. 1. However, the role of 
symmetries in the analogue simulation is subtle as they are ar- 
tificial. The Hamiltonian engineered in the lattice corresponds 
to a certain reference frame (or gauge) only. To consider trans- 
formed Hamiltonians in other frames implies physical modifi- 
cations of the experimental apparatus - for similar discussions 
of gauge symmetries see [27]. In practise, we can choose how 
M transforms. The present bounds on Lorentz violation [28] 
in the accessible energy regime are very tight, and far from 
the current accuracy achievable in optical lattices - the ma- 
jor limitation is due to the lattice correction to the continuum 
theory, hence ~ l/(# lattice points). However, the simula- 
tor is a useful tool as the parameters can be tuned almost at 
will and dynamically generated flavour couplings and further 
strong coupling interactions [29] may be included. 

Helicity projection and mixing with anti-neutrinos. Above, 
we have considered a single Dirac point K = | ( 1 , 1 ) . Since 
the system involves both positive and negative helicities [7], 
we briefly discuss how to achieve helicity projection for com- 
pleteness. We note this issue is closely related to the so called 
valley polarization [30] in the graphene community where 
the degeneracy of the two inequivalent Dirac points related to 
time-reversal symmetry is lifted such that only a single val- 
ley is occupied. A specially designed optical potential barrier 
may be used to filter out one valley following the discussion 
in [31, 32], thus achieving the helicity projection. Alterna- 
tive method would be to engineer a Wilson mass term in the 
optical lattices which decouple doublers thus leaving a single 
Dirac point [33]. In order to consider the mixing between neu- 
trino and anti-neutrino, a Bragg pulse that connects the pos- 
itive and negative helicities may be used as described in the 
(3 + 1) Dirac fermions in optical lattices [34]. An alternative 
approach is to couple a two layered system where the sign 
of the topological charge in one layer is flipped, by sending 



Ty -Ty. 

Engineering exotic particle dispersions. Modifications of 
dispersion relations have attracted a lot of interest in both con- 
densed matter and high energy physics. The Dirac fermions 
in graphene and semi-Dirac fermions with linear dispersion 
in one direction and quadratic dispersion [35] in the other are 
examples of engineered dispersion relations in a condensed 
matter setting. On the other hand, in high energy physics, 
modifications of the energy momentum dispersion relations at 
the Planck scale are suitable, for instance, to reconcile Lorentz 
symmetry and finite resolution of the spacetime points in 
quantum gravity models [12]. As an alternative to the split- 
ting of the multiple degenerate Dirac point, it is possible to 
preserve the Dirac point and to engineer exotic particle dis- 
persion to have any power of momentum. This Dirac-point 
preserving mixing is described by the Hamiltonian Jtf(p) = 
{superdiag{g<7 + },diag{p- <7},subdiag{g<7 - }}, which mixes 
the N Dirac species and gives the spectrum E ~ p N [13]. Our 
setup also allows us to investigate the mixing of Dirac-Weyl 
fermions with high spin beyond the above spin 1/2 case, for 
instance the fate of the flat band with integer spin s after mix- 
ing. Furthermore, disorder can also be introduced to the op- 
tical lattice [36] to simulate spacetime fluctuations at Planck 
scale, thus giving an opportunity to resolve possible conun- 
drums in quantum gravity models by a combination of con- 
cepts, such as modifications of dispersion relations, spacetime 
fluctuations, and Lorentz violations discussed above. 

Conclusions. We have shown how multiple-layered Dirac 
cones can be used to simulate various phenomena from the 
realm of particle physics. The models discussed here deal 
with a quadratic Hamiltonian, but quartic [29] and dynami- 
cal gauge [37] interactions may be included. This, together 
with the possibility of changing the physical parameters in a 
wide range due to the advances in manipulation of cold atoms 
in optical lattices, makes the quantum simulator a useful tool 
to test fascinating ideas from condensed matter to high en- 
ergy physics. Furthermore, such quantum simulations can 
provide a guide to design nano devices, based on novel physi- 
cal effects. We believe the interplay between condensed mat- 
ter, particle physics and quantum simulations could lead us to 
new paradigms, much in the spirit of fundamental theories as 
emerging phenomena [38]. 
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Note added. Once this work was completed, [39] appeared 
proposing an ion trap simulator of NO in 1+1 dimensions. 
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